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Abstract 

If X is a graph with adjacency matrix A, then we define H{t) to be the 
operator exp(/f A). The Schur (or entrywise) product H{t)oH{-t) is a dou- 
bly stochastic matrix and because of work related to quantum computing, 
we are concerned with the average mixing matrix Mx, defined by 

Mx^ lim - / H{t)oH{-t)dt. 

C^cx) C Jo 

In this paper we establish some of the basic properties of this matrix, show- 
ing that it is positive semidefinite and that its entries are always rational. 
We see that in a number of cases its form is surprisingly simple. Thus for 
the path on n vertices it is equal to 

^ {2J+I+T) 



2n + 2 



where T is the permutation matrix that swaps j and n + l- j for each j. If 
X is an odd cycle or, more generally, if X is one of the graphs in a pseudo- 
cyclic association scheme on n vertices with d classes, each of valency m, 
then its average mixing matrix is 

n-m + l m-l 
J+ /. 



(One reason this is interesting is that a graph in a pseudocyclic scheme 
may have trivial automorphism group, and then the mixing matrix is more 
symmetric than the graph itself.) 
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1 Average Mixing 



Let X be a graph with adjacency matrix A. We define a transition matrix Hxit) 
by 

Hxit) := expiitA]. 

For a physicist this matrix determines a continuous quantum walk. We note 
that it is both symmetric and unitary, in particular Hx ( t) = Hx {-t). For relevant 
recent surveys see, e.g., Kendon and Tamon |6|, Godsil |4|. 

The matrix Hit] gives rise to a family of probability densities as follows. Let 
AoB denote the Schur product of two matrices A and B with the same order. 
Thus 

iAoB)u,v^Au,:,Bu,v 

and we will use A°^ to denote Ao A. Then if 

Mxit):^Hxit)oHxi-t) 

we see that Mxit) is a nonnegative real matrix and, since Hxit) is unitary, each 
row and column of Mxit) sums to 1. We use to denote the standard basis 
vector of C^'"^^ indexed by the vertex u, thus we have a family of probability 
densities e^Mxit), and we are concerned with the behavior of these densities. 
In this paper we are interested in the matrix Mx, which we define by 

Mx = - Mxit)dt. 

C Jo 

Following I LI, we call this the average mixing matrix of X. 

To work with Mx and Mx, we use the spectral decomposition of A. This 
allows us to write A as 

r 

where 0,- runs over the distinct eigenvalues of A and is the matrix represent- 
ing orthogonal projection onto the eigenspace belonging to 0, . Then we have 

Hxit) = Y,ewiidrt)Er 

r 

and 

Mx (?) = E ^r' + 2 E COS((0, - Bs) t) Br o 
r r<s 

From this we have: 
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1.1 Lemma. If A = OrE^ is the spectral decomposition ofS = AiX), then 



2 Properties of the Average Mixing Matrix 

In this section we derive some basis properties of the average mixing matrix. 
We know aheady that it is symmetric. By a famous theorem of Schur, the Schur 
product Mo of two positive semidefinite matrices is positive semidefinite and, 
since the sum of positive semidefinite matrices is positive semidefinite, we see 
that Mx is positive semidefinite. If A and B are symmetric matrices of the same 
order, we write A)^ Bio denote that A- Bis positive semidefinite. 
Clearly Mx is a nonnegative matrix. In fact: 

2.1 Lemma. IfX is connected, all entries ofMx are positive. 

Proof. Each Schur square is nonnegative and if [Mx) u,v-0 then {E°^) «, ^ = 0. 
However this implies that (Er) u,v -0 for all r and hence any linear combination 
of the idempotents Ey has wi'-entry zero. Since this implies that iA'^)u,v - for 
all k we conclude that X is not connected. □ 

In m Lemma 16.2] we show that average mixing is never uniform, that is, 
Mx cannot be a scalar multiple of the all-ones matrix /. 

We say that a matrix M is a contraction if x* Mx < x*x for eill complex vec- 
tors X. 

2.2 Lemma. IfMx is the average mixing matrix of the graph X then all eigenval- 
ues of Mx lie in the interval [0, 1] . 

Proof Since H{t) ® H{-t) is unitary and since H{t) o H{-t) is a principal sub- 
matrix of Hit) ® H{- 1), it follows that H{t) o H{- 1) is a contraction. Since it is 
symmetric and real, its eigenvalues must lie in the interval [0, 1], and therefore 
this holds true for Mx too. □ 

Since Mx is the average of doubly stochastic matrices, it is doubly stochastic 
and its largest eigenvalue is 1. We can also see this without appealing to the 
averaging. Notice that 

{{EroErWu = iielEr)o{elEr))l = {elEr,elEr) = e^E^eu = e^ErCu = (£",-)«,« 
and, since = /, it follows that Mxl = 1. 
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2.3 Lemma. The average mixing matrix of a graph is rational. 

Proof. Let (p(.X, x) be the characteristic polynomial of X, and let F be a splitting 
field for (piX, x). We use the fact that an element of F which is fixed by all field 
automorphisms of F must be rational. If a is an automorphism of F, then 

A = A'' = Y,OrK- 

r 

Since 0^ must be an eigenvalue of A and since the spectral decomposition of 
A is unique, it follows that is one of the idempotents in the spectral decom- 
position of A. Therefore the set of idempotents is closed under field automor- 
phisms, and so must the {-6'°^}^. Consequently 

for all a and therefore Mx is rational. □ 

Note that this lemma holds whether we use the Laplacian or the adjacency 
matrbc — all we need is that A be symmetric with integer entries. 

We use L{X) to denote the Laplacian matrix X. If A is the diagonal matrix 
whose i-Xh diagonal entry is the valency of the /-vertex of X, then 

L{X) = A-A. 

The transition matrix of X relative to L{X) is 

Hi,(t):=exp(zf(Z\-A)). 

When X is regular, questions about Hl reduce immediately to questions about 
Hx, but in general there is no simple relation between the two cases. 

2.4 Lemma. IfX is regular then X and its complement X have the same aver- 
age mixing matrix. For any graph X, the average mixing matrix relative to the 
Laplacian ofX is equal to the average mixing matrix relative to the Laplacian of 
X. 

Proof. If X is regular then the idempotents in the spectral decomposition of 
its adjacency matrix are the idempotents in the spectral decomposition of the 
adjacency matrix of X. For any graph X on n vertices 

L{X)^L{Kn)-L(Xy, 

since L{Kn) = nl - J and since LiX) commutes with /, the idempotents in the 
spectral decomposition of its Laplacian are the idempotents in the spectral de- 
composition of the Laplacian of X. □ 
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3 Integrality 



5r,s- 



In investigating the relation between the structure of Mx and the graph X, it 
can be convenient to scale Mx so that it entries are integers. For this we need 
to know a common multiple of the denominators of its entries. 

3.1 Lemma. IfDis the discriminant of the minimal polynomial of A, then Mx 
is an integer matrix. 

Proof Let 6i,...,9m be the distinct eigenvalues of A. Define polynomials £rit) 
by 

erit):=Y\it-es). 

We note that (ridy) = y/'{9r) and 
Now 

Er = -^E-JriA). 

The discriminant D of if/ is equal (up to sign) to 

m 

since the entries of f ^ (^) are algebraic integers we conclude that the entries of 
D^E°^ are algebraic integers and therefore the entries of D^Mx are algebraic 
integers. Since Mx is rational, the lemma follows. □ 

We have no reason to believe this lemma is optimal. If the eigenvalues of X 
are all simple, we can do better. 

3.2 Theorem. Let X be a graph with all eigenvalues simple and let D be the 
discriminant of its characteristic polynomial. Then DMx is an integer matrix. 

Proof We have 

ixi-Ar'^Y.^T^'- 
T x-e,- 

and since (/- tA)~^ is the walk generating function of X, it follows from O 
Corollary 4. L3] that 

(X - Or) [(p{X\ u, x)0(X\ V, X) - (p{X\uv, x)(p{X, x)Y'^ 

iEr)u,u = lim 

<p{X,x) 
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and since 

(l>iX,x) , ^ 
lim^!^-^ =(P'iX,9r) 

we conclude that if 0^ is simple 

^2_ 4>{X\u,er)(p{X\v,er) 

If B is the n X n matrix with ur-entry u,6r) and A is the ny- n diagonal 
matrix with r-th diagonal entry (p'{X, Or), it follows that 

Mx = BA~^b'^. 



Assume n = \ V{X) \ and let 0i, . . . , 0„ be the eigenvalues of X. Let V be the 
nx n Vandermonde matrix with -entry Let (p be the characteristic poly- 
nomial of X. The discriminant of (p is equal to the product of the entries of A, 
we denote it by D. 

Let C be the nx n matrix whose i/r-entry is the coefficient of x^~^ in (p{X\ 
u, x). Then CV = B and 

Mx = CVA~^V^C^. 

Define polynomials ^^(x) by 

£sM = Y\{x-er) 

Tits 

and let L be the nxn matrix with 5 j-entry equal to the coefficient of x-'"^ in 
es{x). Note that 

es{er) = 6r,s(t>'{es) 

and therefore 

LV = A. 

Since 

A = A'^ = V^L^ 

it follows that 
and so 
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(We're almost done.) The entries of L are algebraic integers. As 

detiVV^) = detiVf = D 
and as LD = A, we see that det(L^L) = D. Therefore the entries of 

are algebraic integers. So the entries of DVA~^V^ are algebraic integers and, 
since the entries of C are integers, the entries of DMx are algebraic integers. 
But the entries of DMx are rational and therefore they are all integers. □ 

It is at least plausible that if D is the discriminant of the minimal polynomial 
of X, then DMx is an integer matrix. 



4 Average Mixing on Paths 

We need the following trigonometric identity. 
4.1 Lemma. 

« ^ sm{{n + be) 
2j:cosm= . ' +1. 
r=o sm(i0) 

Proof. Ifq:^ e'^ then 



n 



r=0 r=0 



+ - \ 

q-l q~^-^ 

<7"+^ - 1 q-^-q 

q-l l-q 

q^+^-q-" 



q-l 



+ 1 



qn+ll2 _ q-ll2-n 



2^ +1. □ 
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We will also use the following explicit expression for the idempotents in the 
spectral decomposition of A(P„). This result is standard, but we offer a sketch 
of the proof. 

4.2 Lemma. The idempotents Ei,...,E„ in the spectral decomposition ofPn are 
given by 



iEr)j,k = 



n + l 



■ sin 



sm 



km 
[n + l 



Proof. If A = AiPn) and e„ is the n vector in the standard basis of U", then 



-sin((n + 



' sin()6) ' 




sin(2)6) 




' sin(^) ' 


sin(2^) 




sin()6) + sin(3^) 


= 2cos(j0) 


sin(2/3) 










,sin(n/3), 




, sin((n-l))6) , 




,sin(n/3), 



Soifsin((nlj6)) =Othen 



( sin(^) ^ 
sin (2^) 



[sinin^], 

is an eigenvector for A with eigenvalue 2 cos(/3) . Letting /3 vary over the values 

2nr 



n + l 



r = l,...,n 



we obtain n distinct eigenvalues. Therefore each eigenvalue of P„ is simple and 
the projection onto the eigenspace spanned by z(/3) is 



zil5)Tzil3) 



We can compute the value of the inner product using Lemma l4Jl 

and this yields the stated expression for iJ^ • □ 



4.3 Lemma. If Ei,...,En are the idempotents for Pn, then the average mixing 
matrix ofPn is 

yEroEr = ^—{2J+I+T). 

y 2n + 2 
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Proof. We use Lemma l4[2l 



{E,. o Er)j,k = sin2 ^ sin^ 



km 



which implies that 

{n + iy 



1 ( i2irn\\ I i2krn\ 

■ {Er o Er)j,k = 7 1 - COS 1 - COS 



Now 



f2irn\\( (2krn\\ 



(^]rn\ (2krn\ 1 i2{i + k)rn\ 1 (2U -k)rn\ 

1 - cos - cos + - cos + - cos 

yn + li ^n + li 2 y n + \ > 2 V n + l ' 

We need to sum each of the five terms on the right from 1 to n. From Lemma fO] 
it follows that 



E(2£ni\ i 
cos = - 

,-1 yn + ll 2 



-1 + 



1 + 



1. 



Consequently 



Y,in + lfiEroEr)j,k^< 



3(n + l)/2, j = k; 

3(n + l)/2, j + k^n + 1 

n + l, otherwise 



and this completes the proof. 



□ 



5 Path Laplacians 

Let A denote the diagonal matrix with i-th diagonal entry equal to the valency 
of the i-th vertex of X. Then the Laplacian of X is Z\ - ^. Let D be the n x - 1) 
matrix with 

Di,i = l, Di,i-i^-l 

and all other entries zero. Then D is the incidence matrix of an orientation of 
Pn and 

DD^ = A-AiP„), D'^D = 2I-A{Pn-i). 
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If 01 , ... , 6n-i are the eigenvalues of Pn-i, this shows that the non-zero eigenval- 
ues of Z\ - AiP„) are the numbers 2-dr tor r = I,... ,n- I. We can also use this 
to determine the idempotents in the spectral decomposition of Z\ - A(P„). The 
lemma follows. □ 

5.1 Lemma. IfEi,..., En-i are the idempotents in the spectral decomposition 
ofPn-i, then the idempotents of A - are n~^J and 

1 

DErD\...,r = l,...,n-l. 



2 -Or 
Proof. Since 

DErD^DEsD^ = DEr {21 - AiPn-i)EsD^) 
and since ErEs = if r 5 and (2/ - A{Pn-i)Es = (2 - es]Es, it follows that 

DErD'^DEsD^ = 5r,si2 - er)DErD^. 

Therefore (2 - 0, )~^DErD^ is an idempotent. Next 

{A - AiPn))DE,.D^ = DD^DErD^ = Di2I - A{Pn-i))ErD^ = (2 - er]DErD^ 

and therefore {2-0, DEr represents orthogonal projection onto an eigenspace 
of A{Pn) ■ The lemma follows. □ 

5.2 Lemma. IfEi,..., E„-i are the idempotents in the spectral decomposition 
ofPn-i, then 

^ , T 2 {{2i-l)rn\ ({2k-l)r7t\ 

{2 - 9 r)~\DErD^)j,k = -cos ' \ ' cos ' \, l< ,k<n. 

n \ 2n I \ 2n ) 



Proof. From Lemma FOl we have 

\ I lcrn\ 

1 < i,k<n-\. 



2 

{Er) j,k = - sm 



n 







(krn\ 


sin 




. n ) 




. n j 



Let a = rnin and let a denote the column vector of length n-l where a j 
sm{ia). Then 

DErD'^ ^-Da{Daf 
n 
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and 



Da = 



sin (a) 
sin(2a) - sin(a) 

sin((n- l)a) -sin((n-2)a) 
-sin((n-l)a) j 



= 2sin(a/2) 



f cos(a/2) \ 
cos(3a/2) 

cos((2n-3)a/2) 
Uos((2n-l)a/2); 



where in computing the last entry we have used the fact that na = rn, whence 
sin(na) = and 

-sin(n- l)a = sin(na) - sin(n - l)a. 



Finally for Pn-i we have 

irjt\ ■?trn\ 
2 - 0r = 2 - 2 cos[— J = 4 sin^ J 

5.3 Theorem. The average mixing matrix for the continuous quantum walk us- 
ing the Laplacian matrix of the path Pn is 

Fr = (.2-Grr^DErD'^. 

(i2j-l)rn\ 



Proof. Set 
Then 



1 



= — T 1 1 + cos 
1^ 



2n 
[2}-l)rTi 



cos 



{2k-l)rn 



2n 



n 



1 +COS 



( {2k-Drn 
\ n 



and 



(l.cos(^))(l.cos( 



{2k-l)rn 



n 



1 ({2i + 2k-2)rn\ 1 [i2j-2k)rn\ 
+ - COS + - cos . 

2 V n ) 2 \ n ) 
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From Lemma l4TT] we have 



2 y cos 

r =1 ^ n I 



It is now easy to derive the stated formula for the average mixing matrix. □ 

We note that 21 - L{Pn) can be viewed as the adjacency matrix of a path on 
n vertices with a loop of weight one on each end-vertex. Examples show that 
if we add loops with weight other than or 1, the average mixing matrix is not 
a linear combination of /, / and T. Thus if we add loops of weight 2 to the 
end-vertices of P^, the average mixing matrix is: 





( 


599 


218 


146 


146 


218 


599 ' 






218 


455 


290 


290 


455 


218 


1 




146 


290 


527 


527 


290 


146 


2*9*107 




146 


290 


527 


527 


290 


146 






218 


455 


290 


290 


455 


218 






599 


218 


146 


146 


218 


599 ^ 



(Here the discriminant of the characteristic polynomial is 2^3^107.) 

6 Cycles 

We determine the average mixing matrices for cycles. 

Let P be the permutation matrix corresponding to a cycle of length n and 
let C be a primitive n-lh root of unity. Define matrices Fq, . . . , F„_i by 



-1 + 



-1 + 



= -1 + 



sin(g) 
-cos(^7r)sin(g) 



sin(g) 



= -((-1)^ + 1). 
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(Thus the rows and columns of these matrices are indexed by {0, . . . , n - 1.) Then 
P is a normal matrix and has the spectral decomposition 

n-l 
r=0 

We also note that 

1 

FroFs = -Fr+s 
n 

where the subscripts are viewed as elements of Z„. The adjacency matrbc of the 
cycle on n vertices is P + P^. Define Eq to be Fq and, if < r < n/2, we set 

= Ff + Ffi-j-. 

Further Eq :- Fq and, if n is even then En/z '■- Enii- Then if n is odd, 

'£'o> •••>-£'(«- 1)/2 

are the idempotents in the spectral decomposition of A(C„) and the correspond- 
ing eigenvalues are 

n-l 



er = C+r, r = 0,...,- 



2 ' 

if n is even we have the additional idempotent En/2 with eigenvalue ("^^ = - 1. 
6.1 Lemma. If n is odd then the average mixing matrix of the cycle Cn is 



n-l I 
-^J+-I, 
n 



if n is even it is 



n-2 , 1 



n^ n 



Proof. Assume first that n = 2m+l. Then the average mixing matrix is 

m m 



r=0 r=l 



1 1 ^ 
= -Fo + - ^ (F2r + F_2r + 2Fo) 

n n ~i 

n-l 

= Fo + -Y.Fr 

1 n-l 

= -!+ Fo. 

n n 
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Now suppose n = 2m. Then the average mixing matrix is 

m-1 m-1 
^m + E Ef = + F°^+J: [Fr + F^rf' 
r=0 r=l 

1 1 1 

= -i^o + -Fo + - ^ (F2r + F_2r + 2Fo) 
n n n fT[ 

m-1 

= Fo + - ^ iF2r+F-2r) 

n-2 2"^^ 

= Fo + - E F2r. 

n ^ r^o 



Since 



we see that P"^ has spectral decomposition 

n-l 

p'^=Y^i-irFs 

5=0 

and consequently 

1 m-1 

Ul + prn^=Y^F2r. 

^ r=0 

Our stated formula for the average mixing matrix follows. □ 

In general the average mixing matrix for a circulant can be more complex 
in structure than the average mixing matrix of a cycle. We will see in the next 
section that graphs in pseudocyclic schemes provide the right generalization of 
cycles or, at least, of odd cycles. 

It is easy for different circulants of order n to have the same spectral idempo- 
tents, any two such circulants necessarily have the same average mixing matrix. 



7 Pseudocyclic Schemes 

An association scheme ^ with d classes on v vertices is a set of 01 -matrices 
Aq, ...,A(i such that 

(a) Ao = IandZiAi=J. 
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(b) For each / we have aT e ^. 

(c) AiAj = AjAi. 

(d) For all i and j the product AjAj lies in the real span of s^. 

Convenient references are [2l|3]. We often view the matrices At as the adja- 
cency matrices of directed graphs on v vertices, we call these the graphs of 
the scheme and note that the axioms are often stated in terms of the graphs 
rather than their adjacency matrices. An association scheme is symmetric if 
each matrix in it is symmetric. These axioms imply that the span of .sa?' is a com- 
mutative matrix algebra that is closed under the Schur product. It is called the 
Bose-Mesner algebra of the scheme. The matrices in .$2/ form a basis consisting 
of Schur idempotents. There is also a basis of matrix idempotents EQ,...,Ed, 
which satisfy 

(a) Eo = v-^JandZEj = I. 

(b) For each j there is an index / such that Ej = Ejf. 

(c) For all i and j the Schur product Ei o Ej lies in the Bose-Mesner algebra. 

The matrix idempotents can be viewed as providing (an analog of) the spectral 
decomposition. 

Since the matrices Ai commute it follows that Ai and / commute and hence 
there are constants Vi such that AJ = viJ, these are the valencies of the scheme. 
The column space of Ej is an eigenspace for all matrices in the Bose-Mesner 
algebra, its dimension is denoted by rrij and it is referred to as a multiplicity 
of the scheme. Note that vq - rriQ - 1. K scheme is pseudocyclic if mi, ... , ma 
are equal (in which case their common value is [v - I)/ d). If a scheme is pseu- 
docyclic then vi,..., Vd are necessarily all equal to {v - \)l d. For details see 
Brouwer, Cohen and Neumaier|2, §2.2B]. 

We note one class of examples, the cyclotomic schemes. Assume q is a 
prime power and d divides q - I. Let F be the finite field of order q and let S 
be the subgroup of the multiplicative group of F generated by the non-zero d- 
th powers. Thus |S| = (^7 - Let Si,...,Srf denote the cosets of Sin F*, with 
Si = S. Now we define the matrices of an association scheme with d classes and 
with vertex set F by setting Aq- I and 

(A)x,y = l, {i = l,...,d) 
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if and only if y - x is in Si. These matrices form the cyclotomic scheme with d 
classes on F. It is symmetric if and only if - 1 £ S. The directed graphs Xi , . . . , 
such that Ai - A{Xi) are all isomorphic. The most well known case is when 
d = 2 and ^7 = 1 modulo four, in which case the two graphs we have constructed 
are the Paley graphs. 

We note that there are pseudocyclic schemes which are not cyclotomic, and 
that there are pseudocyclic schemes with two classes where the two graphs are 
asymmetric, that is, their only automorphism is the identity. 



8 Average Mixing on Pseudocyclic Graphs 

8.1 Theorem. Suppose X is a graph in a d -class pseudocyclic scheme on n ver- 
tices consisting of graphs of valency m = {n-\)ld. Then the average mixing 
matrixofX is 

n - m + l m-l 
5 /+ /. 

Proof Koppinen LZI proved that for an association scheme with d classes we 
have 

d d Y 

If the scheme is symmetric then Ai o At = At and, since Mo N is a submatrix of 
N, Koppinen's identity yields that 

^1 ^1 
tonvi p^mj 

For any scheme, we have mo = vq = I and for a pseudocyclic scheme 

mi = Vi = ——, i^l,...,d. 
d 

As = J and Y.] Ej = I,ifm = {n-l]/dwe find that 

m ' Jtt- p ' 



and hence 



„ (I l\ m-l n^-nm m-l 

^ /+ 1 = ^ /+ 1. n 



-LEf 



nm n^ } nm n-'m nm 
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A graph in a pseudocyclic scheme with two classes is known as a conference 
graph, and there are examples whose automorphism group is trivial. 

9 Cospectral Vertices 

If M e V[X) then (A*^) 

u,u is equal to the number of closed walks on X of length 
k that start (and finish) at u. The following result summarizes some standard 
facts. (The equivalence of (a) and (c) follows from Equation (1) on page 52 of 
m and (a) and (b) from Equation (2) on page 9 of the same source.) 

9.1 Theorem. Let X be a graph with vertices u and v, and let Ei,...,Em be the 
idempotents in the spectral decomposition of AiX). Then the following are 
equivalent: 

(a) iA'']u,u = iA'')^^^forallk>0. 

(b) \\Ereuf^\\Ere,f. 

(c) The graphs X\ u and X\ v are cospectral. □ 
Note that 

We will say that vertices u and v oiX are cospectral if any of the conditions in 
this theorem hold. 

We recall that the graph X admits perfect state transfer from u to at time 
T if there is a scalar j such that 

where lyl = 1 since 11^^11 = \\Hx{T:)eu\\- Hxit) is a polynomial in A we have 

ErHxit) = Hx{t)Er = exp[iert)Er 

and therefore 

yErei; = Hxir)Ereu = exp(/0,-T)£'re,;. 

Since ErBi, and ErCu are real vectors and since lyl = |exp(i0rf)l = 1> we con- 
clude that 
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(This fact is well known and we include the proof for convenience only.) 
Clearly if ErCu = ±£^6^ then 

WEreuf =\\Ereyf- 

and so u and v are cospectral. We will say that u and v are strongly cospectral 
if Ereu = +Erev. 

9.2 Lemma. If u and v are vertices in X and the eigenvalues of X are simple, 
then u and v are strongly cospectral if and only if they are cospectral. 

Proof We only need to consider the case where u and v are cospectral. The 
vectors -Er^w and EyBi, are both eigenvectors of X with eigenvalue Or- Since Oy 
is simple this means that that E^eu is a scalar multiple of E^eu- Since these two 
vectors have the same length, the lemma follows. □ 

9.3 Theorem. Let Mx be the average mixing matrix of the graph X. Then ver- 
tices u and V are strongly cospectral if and onlyifMxieu - e^] = 0. 

Proof Suppose iV )>= and Niei - 62) = 0. We may assume that the leading 2x2 
submatrix of N is 

(a b\ 
[b d] 

and therefore 

= {ei- 62)^ N {61-62) = a + d-2b. 
Hence b = ia + d)/2. Since AT ^ we have ad-b^ >0 and thus 

< 4ad - 4b^ = 4ad -ia + df - -[a- df, 

whence a = d = b. 

IfMx(e„-eJ = Othen 

r 

and as each summand E'f' in Mx is positive semidefinite, we have 

Efi6u-ei;] =0 

for all r. Therefore, for all r, 

aEr)u,uf = iiEr)u,vf = iiEr)v,vf 

Since 

iEr)u.v = {Er6u,Erev) 

it follows by Cauchy-Schwarz that E^eu = i^'r^y □ 
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A graph is walk regular if all its vertices are cospectral. There are many 
classes of examples, and it is not necessary that such graphs be vertex transi- 
tive. We note the following though: 

9.4 Lemma. If all vertices in X are strongly cospectral, then X is Ki or 

Proof. If all vertices are strongly cospectral, then all rows of Mx are equal, 
and therefore it is a scalar multiple of /. By |4, Lemma 16.2] this implies that 

\ViX)\<2. □ 

One consequence of the results here is that, if the rows of the average mix- 
ing matrix of X are distinct, there is no perfect state transfer on X. Kay |5, Sec- 
tion D] has shown that if we have perfect state transfer in X from utov and 
from u to w, then v-w. 



1 Discrete Walks 

If [/ is a unitary matrix then its powers describe a discrete quantum walk. The 
corresponding average mixing matrix is 

Af-l 

lim y VoU'". 

Since unitary matrices are normal, they have a spectral decomposition which 
we can write as 

r 

where 0^ is real and Ey Hermitian for all r. The unitary matrices considered in 
quantum computing are constructed from an underlying regular graph. 

10.1 Theorem. lfE\,...,Em are the idempotents in the spectral decomposition 
of the unitary matrix U, the average mixing matrix UofU is Y.r Ef'- If U is 
rational then U is rational. 

Proof We have 

[/« o [/-« = J^Ef + ^ o Es. 

r r^s 

Therefore 



1 iV-l 1 lN-\ \ 



s 
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and since 

N-l „ni{er-6s] _ 1 
y ni[Br-e,] ^ f ^_ 

«=o ^ - 1 

The right side is bounded in absolute value by 

2 



_ i| 



which is independent of N. Hence 

Y N-l 



o2 



as N tends to infinity. □ 
Since 

/ - [/" o [/-" = ^ (1 - o = 2 ^ (1 - cos(n(0, - o 

we see that I - oU~^ is positive semidefinite and hence 7 - LA is positive 
semidefinite. 

If U is based on a regular graph then U is rational. 



1 1 Questions 

We list three questions raised by our work. 

1 . Is it true that if D is the discriminant of the minimal polynomial of X, then 
DMx is an integer matrix? 

2. Is there a useful algorithm for computing Mx that works over the ratio - 
nals? 

3. Which graph have the property that their average mixing matrix is a linear 
combination of / and /? (Jamie Smith, in unpublished work, has found a 
strongly regular graph with this property that is not pseudocyclic.) 
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